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Light- front Hamiltonian methods are being developed to attack bound-state problems in QCD. In 
■ this paper we advance the state of the art for these methods by computing the well-known Lamb shift in 
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hydrogen starting from first principles of QED. There are obvious but significant qualitative differences 
between QED and QCD. In this paper, we discuss the similarities that may survive in a non-perturbative 
QCD calculation in the context of a precision non-perturbative QED calculation. Central to the discussion 
are how a constituent picture arises in a gauge field theory, how bound-state energy scales emerge to guide 



^ . the renormalization procedure, and how rotational invariance emerges in a light-front calculation. 

m 

^ ; PACS number(s): ll.lO.Ef, 12.20.Ds 

I. INTRODUCTION 

Why is the calculation of the Lamb shift in hydrogen, which at the level of detail found in this 
paper was largely completed by Bethe in 1947 [Q, of any real interest today? While completing such a 

X" . . . 

' calculation using new techniques may be very interesting for formal and academic reasons, our primary 

motivation is to lay groundwork for precision bound-state calculations in QCD. The Lamb shift provides 
an excellent pedagogical tool for illustrating light-front Hamiltonian techniques, which are not widely 
known; but more importantly it presents three of the central dynamical and computational problems 
that we must face to make these techniques useful for solving QCD: How does a constituent picture 
emerge in a gauge field theory? How do bound-state energy scales emerge non-perturbatively? How does 
rotational symmetry emerge in a non-perturbative light-front calculation? 

These questions can be answered in detail in QED. The answers clearly change in QCD, and we point 
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out several places where this is clear, but we hope that much of the computational framework successfully 
employed in QED will survive. 

In order to formulate these questions in a more precise fashion, we first outline the general compu- 
tational strategy we employ. First, we use the renormalization group to produce a regulated effective 
Hamiltonian H\ , where A is a cutoff and renormalization is required to remove cutoff dependence from all 
physical quantities. At this point we have a regulated Hamiltonian that contains all interactions found 
in the canonical Hamiltonian, a finite number of new relevant and marginal operators (each of which 
contains a function of longitudinal momenta because longitudinal locahty is not maintained in light-front 
field theory), and an infinite number of irrelevant operators as would occur in any cutoff theory. This com- 
plicated Hamiltonian cannot be directly diagonalized, and since we want to solve bound-state problems 
we cannot solve it using perturbation theory. The second step is to approximate the full Hamiltonian, 
using 

H^ = Ho + {Hx-Ho) = Ho + V , (1) 

where Ho is an approximation that can be solved non-perturbatively and V is treated in bound-state 
perturbation theory (BSPT). The test of Ho is whether BSPT converges or not. 

We can now reformulate the questions above. Is there a scale A at which Ho does not require particle 
emission and absorption? What are the few-body interactions in Ho that generate the correct non- 
perturbative bound-state energy scales? Is there a few-body realization of rotational invariance; and if 
not, how does rotational symmetry emerge in BSPT? We should emphasize that for our purposes we 
are primarily interested in answering these questions for low-lying bound-states, and refinements may be 
essential to discuss highly excited states or bound-state scattering. 

It is essential that A, which governs the degree to which states are resolved, be adjusted to obtain a 
constituent approximation. If A is kept large with respect to all mass scales in the problem, arbitrarily 
large numbers of constituents are required in the states because constituent substructure is resolved. 
A constituent picture can emerge if high free-energy states couple perturbatively to the low free-energy 
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states that dominate the low-lying bound-states. In this case the cutoff can be lowered until it approaches 
the non-perturbative bound-state energy scale and perturbative renormalization may be employed to 
approximate the effective Hamiltonian. In QED we note that the range into which the cutoff must be 
lowered is 

ma^ <C A < ma , (2) 

where A = A — rrip — rUe as will be explained later, and m is the reduced mass of hydrogen. If the cutoff 
is lowered to this range, hydrogen bound-states are well approximated using proton-electron states and 
including photons and pairs perturbatively. 

It is an oversimplification to say the constituent picture emerges because the QED coupling constant 
is very small. Photons are massless, and regardless of how small a is, one must in principle use nearly de- 
generate bound-state perturbation theory that includes extremely low energy photons non-perturbatively. 
This is not required in practice, because the Coulomb interaction which sets the important energy scales 
for the problem produces neutral bound-states from which long wavelength photons effectively decouple. 
Because of this, even though arbitrarily small energy denominators are encountered in BSPT due to 
mixing of electron-proton bound-states and states including extra photons, BSPT can converge because 
emission and absorption matrix elements vanish sufficiently rapidly. 

The well-known answer to the second question above is the two-body Coulomb interaction sets the 
non-perturbative energy and momentum scales appropriate for QED. We have already used the results of 
the Bohr scaling analysis that reveals the bound-state momenta scale as p ~ ma and the energy scales as 
E ^ ma^. As a result the dominant photon momenta are also of order ma, and the corresponding photon 
energies are of order ma. This is what makes it possible to use renormalization to replace photons with 
effective interactions. The dominant photon energy scale is much greater than the bound-state energy 
scale, so that A can be perturbatively lowered into the window in Eq. (^) and photons are not required 
in the state to leading order. A similar analysis in QCD will reveal qualitatively different results. If a 
constituent picture emerges, the dominant interaction will be confining and the dominant giuon energy 
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scale will be directly affected by confinement. A confining interaction automatically generates a mass 
gap for gluon production. 

Finally we discuss rotational invariance in a light-front approach. In light-front field theory, boost in- 
variance is kinematic, but rotations about transverse axes involve interactions. Thus rotational invariance 
is not manifest and all cutoffs violate rotational invariance in light-front field theories. In QED it is easy 
to see how counterterms in H\ arise during renormalization that repair this symmetry perturbatively; 
however, the issue of non-perturbative rotational symmetry is potentially much more complicated. We 
first discuss leading order BSPT and then turn to higher orders. 

To leading order in a constituent picture we require a few-body realization of rotational symmetry. 
This is simple in non-relativistic systems, because Galilean rotations and boosts are both kinematic. In 
QED the constitutucnt momenta in all low-lying bound-states are small, so a non-relativistic reduction 
can be used to derive Tio- Therefore to leading order in QED we can employ a non-relativistic realization 
of rotational invariance. This type of approach can be tried in QCD, but it is not essential that it work 
because alternative few-body realizations of the full set of Lorentz symmetries exist. 

At higher orders in BSPT rotational invariance will not be maintained unless corrections are regrouped. 
We have computed hyperfine structure and shown that terms from first-order and second-order BSPT 
are required to obtain angular momentum multiplets The guiding principle in this and all higher 
order calculations is to expand not in powers of V, but in powers of a and log(Q!). Tio should provide the 
leading term in this expansion for BSPT to be well-behaved, and subsequent terms should emerge from 
finite orders of BSPT after appropriate regrouping. Powers of a appear through explicit dependence 
of interactions on a, and through the dependence of leading order eigenvalues and eigenstates on a 
introduced by interactions in Tio- This second source of dependence can be estimated using the fact that 
momenta scale as ma in the bound-state wave functions. Of more interest for this paper is the appearance 
of log(a), which is signaled by a divergence in unregulated bound-state perturbation theory. As has long 
been appreciated, such logarithms appear when the number of scales contributing to a correction diverges. 

The existence of a small parameter simplifies the non-perturbative calculation of bound-state observ- 
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ables considerably, and it has been suggested that a similar expansion be employed to guide light-front 
QCD calculations even if it requires the introduction of masses that violate rotational invariance away 
from the critical value of the coupling Q . We do not detail this proposal, but a thorough understanding 
of such expansions in QED is almost certainly necessary before one has any hope of using this approach 
for QCD. 

We proceed with a description of our Lamb shift calculation. In hydrogen there is a small amplitude 
for a bound electron to emit and re-absorb a photon, which leads to a small shift in the binding energy. 
This is the dominant source of the Lamb shift, and the only part of this shift we compute in this paper. 
This requires electron self-energy renormalization, but removal of all the bare cutoff. A, dependence 
requires a complete 4th order calculation, which is beyond the scope of this paper. We work with a finite 
bare cutoff: A — ?n-\/2, and show that our results are independent of the effective cutoff, A. 

The energy scale for the electron binding energy is mo? , while the scale for photons that couple to 
the bound-states is ma. This energy gap makes the theory amenable to the use of effective Hamiltonian 
techniques. For simplicity, we use a Bloch transformation Q in this paper to remove the high energy 
scale (i.e., ma) from the states, and an effective Hamiltonian is derived which acts in the low energy 
space alone. This effective Hamiltonian is treated in BSPT, as outlined above. The difference between 
the 25*1 and the 2Pi energy levels, which are degenerate to lowest order, is calculated. 

We divide the calculation into two parts, low and high energy intermediate photon contributions. The 
low energy photons satisfy |k| < A, while the high energy intermediate photons satisfy A < |k| < m. A 
is the effective cutoff for the theory, which is chosen to lie in the range given in Eq. (^). This choice lies 
between the two dominant energy scales in the problem and allows us to avoid near degeneracy problems. 
When an actual number is required we use 

A = a\fa TO ~ 6 X 10^'' to . (3) 

Note that the spectrum of the exact effective Hamiltonian is independent of A, but our approximations 
introduce A-dependence. The range for A is chosen so that the errors appear at a higher order in a than 
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we compute. 

One further introductory comment, the high photon energy (A < |k| < to) part of the shift is further 
divided into two regions, A < |k| < 6 and fo < |k| < m, where b is an arbitrary parameter chosen in the 
range ma <^ b <^ m. This simphfies the calculation with appropriate approximations being used in the 
respective regions. The result must obviously be independent of this arbitrary division point 6, and is, 
unless "non-matching" approximations are used in the respective regions. 

We now outline the paper. In §11 we discuss the theoretical framework of this light-front Hamiltonian 
approach, and in §111 we proceed to discuss the origin of the Coulomb interaction in this framework. §IV 
contains the heart of the Lamb shift calculation. In the final section, §V, we summarize and discuss our 
results. 

II. THEORETICAL FRAMEWORK 

In this paper, the proton will be treated as a point particle. The Lagrangian for the electron, proton, 
and photon system is (e > 0) 

jC = —Ff.^F'"' +lp^{i ^ + e4-me)i^e+i^pii ^-e4-mp)i;p . (4) 

The reduced mass of the system is defined in the standard way 



TO 



m,,mp 



TOe (1 - TOe/TOp + ^(I/to^)) . (5) 



TOe + nip 

Note that in this paper we take the limit rup — > oo because we are only interested in the dominant part 
of the Lamb shift. The Lagrangian leads to the following canonical Hamiltonian in the light-cone gauge, 
A+ =0, 



H 



j d^x^dx- n , (6) 
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+ {icy'd' - ea'A' - irup) [{ia'd' - ea'A' + irup) ^p 



-Lj+^J+ + J+^A^ . (7) 

Note that i = 1,2 only; J+ = 2e (C^^p - CiCe) , and cr* are the standard SU{2) PauU matrices. The 
dynamical fields are A\ and ^p, the transverse photon and two-component electron and proton fields 
respectively. For the relation between "0 and ^ and a comprehensive summary of our light-front conven- 
tions see Appendix A. 

The free Hamiltonian is 

h = i?|(e=o) = / E ibl{p)bs{p) P—±!^ + Bl{p)BM P±_±II^ + al[p)aM ^) , (8) 

plus the anti-fermions. The notation for our free spectrum is h\i) = ei\i) with J^i ~ where the 
sum over i implies a sum over all Fock sectors, momenta, and spin. Next, we normal-order all interactions 
and neglect zero modes. The canonical interactions from Eq. that we use in this paper are 

Vie = J d^x^dx^ Vie , "^ip = J d^x^dx^ Vip , V2 — J d^x^dx^ V2 , (9) 

with 

Vie = eCl^M'^ [{la'd' + ime) C^] + e^l (la'd' - ime) [aM'^e] - 'i<%^A' , (10) 

Vip = -e^taM'-i^ [{za'd' -f im^) ^ - {za'd' - imp) [a^A^^ + 2ei%^A^ , (11) 

V2 = -ij+^^J+ . (12) 

These are the photon emission and absorption by the electron, photon emission and absorption by the 
proton, and instantaneous photon interactions respectively. 

Given the canonical Hamiltonian, H , we cut off the theory by requiring the free energies of all states 
to satisfy 

e^ < -^pT— ' (13) 

where A is the bare cutoff, and V = (V'^,V^) is the total momentum of the hydrogen state. Then, with 
a Bloch transformation we remove the states with free energies satisfying 

—pT- ^^^^ r+ ' 



where A is the effective cutoff. The resuh is an effective Hamihonian, acting in the low energy 
(si < ^ p+ ) space alone. We do not discuss the derivation of H\ any further, but instead refer the 
interested reader to Appendix B. 

Given we then make the following division 

Hx=no + [Hx -Ho) = Ho + V , (15) 

where Tio is an approximation that can be solved non-perturbatively (for this QED calculation) and 
V is treated in BSPT. The test of Ho is whether BSPT converges or not and closely related: is the 
A-dependence of the spectrum weakened by higher orders of BSPT? 

III. LOWEST ORDER SCHRODINGER EQUATION 

The primary assumption we make in this QED bound-state calculation is that the Coulomb interaction 
dominates all other physics. In this work we will treat the Coulomb interaction between the electron and 
proton to all orders in all Fock sectors. After this assumption, the kinematic length scale of our system 
is fixed, 

1 1 137 

p ma m 

which then fixes our dynamical time and length scale, 

^ 1 1 1372 

p2/(2m) ma'^ m 

as is well known, dynamical changes occur very slowly in this system. Note that in this QED calculation 
we will treat photons as free since they carry no charge and interact very weakly at low energies. This 
of course changes drastically for QCD since gluons do carry color charge and interact strongly at low 
energies. After choosing Ho, the a-scaling of our BSPT is fixed, and the spectrum is then calculated to 
some desired order in a and log a. 

In the Coulomb gauge, the Coulomb interaction appears directly in the canonical Hamiltonian, which 
of course is not true in the light-cone gauge.0 In the light-cone gauge, the Coulomb interaction arises 

^However, a confining potential does appear directly in the canonical Hamiltonian in the light-cone gauge, which is a 
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from a combination of two types of interactions in our effective Hamiltonian, instantaneous photon 
exchange and the two time orderings of dynamical photon exchange. Graphically this is shown in Fig. 1. 
These interactions arise from first and second-order effective interactions respectively. See Eq. (|l35| ) of 
Appendix B for the form of the effective Hamiltonian, Hx. 

The time-independent Schrodinger equation in light-front coordinates that the sum of the three time- 
ordered diagrams in Fig. 1 satisfies is ^ 

X Vc ^NixKSeSp) . (16) 

Al^ is the mass squared eigenvalue of the state (j>N, where "N" labels all the quantum numbers of this 
state. The tildes will be notationally convenient below. We have introduced the following Jacobi variables 



Pe = ixr+,K + xr^), (17) 

p', = {x'V+ , n' + x'V^) , (18) 
where pe and p'^ are the initial and final electron three-momentum respectively, and 

Pe+Pp=p'e+Pp^V^ir+,V^) (19) 

is the total momentum of the hydrogen state. Note that k is a two-vector. The norm is defined by 

J (i^K/(27r)^y dx/{4:TT) (f)'^(xKSeSp)4>N' (xKSeSp) = Snn'- (20) 

Vc is the sum of the interactions given by the three diagrams in Fig. 1, and will not be written in all its 
gory detail.^ The leading order term of Vc in a non-relativistic expansion is defined as Vc and is written 

convenient starting point for QCDs+i jsj. 

^For a derivation of Eq. ( |l6| ) from the Schrodinger equation in Fock space see Eqs. (81) to (83) in §111. B.l of Ref. [Q. 
^The interested reader should consult Eqs. (70) and (71) and the discussion below in §111. A. 2 of Ref. [Q; these equations 

are for the equal mass case, but are readily generalized to the unequal mass case; also note that in this reference we used 

a similarity transformation instead of a Bloch transformation; the Bloch transformation was chosen for the current paper 

because of its simplicity. 
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below. 

The non-relativistic expansion is defined in the foUowing way. A coordinate change which takes the 
range of longitudinal momentum fraction, x £ [0, 1] to G [—00, 00] is defined: 



This step can be taken for relativistic kinematics, but there may be no advantage. Then, the non- 
relativistic expansion is an expansion in |p|/m; i.e., we assume 

m » IpI , (22) 

where we have defined a new three-vector in terms of our transverse Jacobi variable, k, and our new 
longitudinal momentum variable, Kz, which replaces our longitudinal momentum fraction, x, 

P = (23) 



Note that the free mass squared in the Schrodinger equation, Eq. (y^, after this coordinate change, 
becomes 



p 



l-x 



(m. + m,f + 2(m. + m,) (i^ - P?:!,!;,^^' + O f ) , (24) 



which is invariant under rotations in the space of vectors p — but not invariant under pz boosts. Here we 
begin to see longitudinal boost invariance being replaced by an expanded kinematic rotational invariance 
in the theory, m is the reduced mass given in Eq. (^. 

Now note that the leading order term of Vc in an expansion in |p|/m is contained in 



Vc ~ (me + mpf J- + -^^Oh ) 6s^s'Js,s'^eL , (25) 



where 



q = P'-P (26) 
h = (^A2 - (Vm2 + p2 + + p^yy |^A2 - (^^m^ + p'^ + + p,2 ^ ^ ^ (27) 
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(me + nip) +2{me + mp)- — r - 



X 6l(^A^- (me + mp)^ + 2(me + mp)2j^jj . (28) 

Note that Qh and da are the constraints that arise from the Bloch transformation. 

It is convenient to define new cutoffs which subtract off the total free constituent masses of the state 

A = A - (me + mp) , (29) 
A = A - (me + mp) . (30) 

In the hmit — > oo we require A and A to be held fixed. Note that this implies 

A^ - (me + nipf ^ - _^ A^ (mp^oo) - 

2(me + mp) 2(me + m^) ' 

A^ - (me + m.pf ^ h? (mp-^oo) ^ 

2(me + rup) 2(me + m^) 

In terms of these new cutoffs, 9l and Or above become 



To see the Coulomb interaction arising from the \ep) sector alone, we make the following requirements 
(which are motivated from the previous two equations) 

«: A «: Ipl and A > Ipl , (35) 

m 

also demanded for |p'| of course. These constraints will be maintained consistently in this paper. Given 
these restrictions we have 

Ol ^ I, (36) 
9h ^ I. (37) 

Vr becomes 



Vc ~ (38) 
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where 

^2^ 4eV'\, , 

Vc = (me + nip) I --2- + I Ss^s'Js^s'^ 

= -{me+mp f (^^^ Ss^s'Js^s'^ ■ (39) 

To finish showing how the Coulomb interaction arises in a light-front Hamiltonian approach, we need 
to know the Jacobian of the coordinate transformation of Eq. ( pT] ) , 



"'^^ ^/p^ + m'^Jp'^ + ml ( Vp^ + + \/p^ + "tp 
1 / /I 1 \ (p2 + 2k2) ^|pl3 



It is also convenient to redefine the norm 



1 + ^^ — -— - o^. ^/^ • (40) 



i^TVAT' = / d^K/(27r)^ / (i2;/(47r) 0^(a;KSeSp)0Ar'(a:KSeSp) 

= ^ d^P -^(p)/(167r^) ^^^r(pSeSp)0Ar'(pSeSp) 

= ^ J (fp (j)l[{pSeSp)(j}N'{pSeSp) ■ (41) 

In this last line the tildes are removed from the wave functions by defining 



'N(PSeSp) = \/ Y^'i'NipSeSp) . (42) 

Putting it all together, the leading order expression for Eq. ( p^ in an expansion in |p|/m given the 
restrictions of Eq. ( |35| ) is 

[-fin + ^yNip's',s'p) = £2 J J^^yMps'Jp) , (43) 

which we see is the non-relativistic Schrodinger equation of hydrogen, m is the reduced mass and — /3„ is 
the binding energy defined by 

Pn = ^\'-(^2''7^^ ■ (44) 
Z(TOe + nip) 

The well known bound spectrum is 
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Pn = 5- , (45) 



where Ryd = ma /2 of course. Note that Eq. (|4^) fixes the a-scahng of |p|: 

IpI - ma . (46) 
Thus we see that the restrictions of Eq. ( ^5| ) become 

ma^ <C A <C ma and A ^ ma , (47) 

which is consistent with Eq. ^ as advertised earUer. 

IV. LAMB SHIFT CALCULATION 

Given our lowest order spectrum, we proceed with BSPT. As advertised, this wiU be divided into 
low and high intermediate photon energy calculations. Before proceeding with these respective calcula- 
tions, we discuss whether Coulomb exchange can be treated perturbatively or non-perturbatively in the 
respective regions. 

For the low energy intermediate photon, the Coulomb interaction between the intermediate electron 
and proton must be treated non-perturbatively, whereas this interaction can be treated perturbatively for 
the high energy intermediate photon contribution. This is seen by noting that each additional Coulomb 
exchange contributes a Coulomb matrix element and an energy denominator which is dominated by the 
larger photon energy scale. Thus each additional Coulomb exchange contributes 

ma^ 



For the low energy photon contribution, in principle |k|mi„ = 0, and each additional Coulomb exchange 
can contribute 0(1), and therefore must be treated non-perturbatively. Of course, when the Coulomb 
interaction is treated non-perturbatively, low-energy intermediate protons and electrons form bound states 
from which long-wavelength photons decouple. This non-perturbative effect leads to |k|mi„ ~ 16.64 Ryd; 
see Eq. ( |121| ) below. For the high energy photon contribution, |k|min = A and from Eq. (^ each additional 
Coulomb exchange thus contributes at most 
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^ = V^~8.5 X 10-2, (49) 
A 

and can therefore be treated perturbatively. 

A. Low energy contribution 

The low energy shift arises from two sources which are shown in Fig. 2. The first term comes from the 
fow energy photon emission part of the effective Hamihonian, {a\vic\b), treated in second-order BSPT. 
Recah Eqs. (^) and (|l^) for the form of wie-Q The second term is the resuh of renormalizing the one 
foop electron self-energy: a counterterm is added to the second-order self-energy effective interaction in 
{a\H\\b), which results in a finite (except for infrared divergences) shift to the electron self-energy. This 
is shown in Fig. 3. The counterterm is fixed by requiring the electron self-energy to evolve coherently 
with the cutoff. The details of defining this counterterm, for the equal mass case, can be found in §111. A. 1 
of Ref . 1^ . A discussion of the physical ideas behind coupling coherence can be found in Ref. |]. For 
further comments on coupling coherence, see the paragraph containing Eq. (39) of Ref. 

Before proceeding with the calculation, we define the binding energy of hydrogen, —B^, in terms of 
the mass-squared, Af^, 

= (me + mp + BnY • (50) 

Assuming Bn is finite as rrip — > oo we have 

^N-i-^^^+^^pf ^ B. + -^^'-^-'b,. (51) 
2(me -|- mpj 2(me + nip) 

Recalling Eq. (^), which is the definition of the zeroth order binding, — in terms of the zeroth order 
mass-squared, and also defining the mass squared corrections, 5M^, by 

Ml, = Ml + SM%, (52) 



* Note that the term where the proton emits and subsequently absorbs a photon is down by two powers of the proton 
mass with respect to the term where the electron emits and absorbs a photon. This result is subtle though, because it is 
true only after the light-front infrared divergences have canceled between two diagrams analogous to the ones in Fig. 2. 
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combined with Eq. (|5l|), gives 



B. = + + (53) 



2 (me + nip) 2mp 
Defining the binding corrections, —SB^^, by 



combined with Eq. (d3), gives 



Bn = Pu + SBn, (54) 



a useful formula to be used below. This formula is useful because SMf^ is calculated below, but SBjsf is 
the quantity that is measured. 

The low energy calculation proceeds as follows. The first term of Fig. 2 is a second-order BSPT shift 
which contributes the following to the mass squared eigenvalue: 

^^ii = E/E^ n;^.n/.,^2 ' (56) 



A" 



k 



DENi{Volf 

where k and are the photon's three-momentum and spin respectively, V = (7'^,7'^) is the total 
momentum of the hydrogen state , and vie is the photon emission interaction given in Eq. (^. Oli 
restricts the energies of the initial, intermediate and final states to be below the effective cutoff ^4. — . 
The explicit form of these restrictions is discussed below. Continuing the description of Eq. ( |5^ ) , 

r _ f <ek^dk+e{k+) _ r d^k 

~ J le^r^F " J (27r)3(2|k|) ' ^ > 

The last step comes from recalling that for a photon fc+ = fc" + fc'^ = |k| + k^. The denominator and 
volume factors are 

Vol = J d^x^dx- = 167r3 6^ {V ~V) , (58) 

The two-body states are 
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\^n{V)) = / ^/pfp+16n^^^r-pe-Pp)^N{PePp)bliPe)Blipp)\0) , (60) 

\^N'{V-k)) = / ^k+k+WTT^S^ (p-k~ki- k2) 4>N'{hk2)bl,^{h)Bl,^{k2m , (61) 



fclfe2 

where (j>N are solutions to Eq. (^3|), the non-relativistic Schrodinger equation of hydrogen, and c^jv is 
related to cj^N by Eq. (||). 

Straightforward algebra leads to 

6Ml, = I I H-P^-Pt) I Oh^-p'^-) f (ptp'e^ktkt) 

X (I67r3(53(fc + A:3-Pe)) (l6^3(53(fc + fci - p^)) 0^ (p^, P - p'J 

X0jV,(fcl,P-fc-fci)0^, (fc3,7'-fc-fc3)0iv(Pe,7'-pe)^i^ , (62) 

where N and A'^' are shorthands for (n, m;) and (n', Z', mj) respectively, the usual principal and angular 

momentum quantum numbers of non-relativistic hydrogen. The "c" on the sum emphasizes the fact that 

the continuum states must be included also. See Eq. (|5|) for DEN\. N\ is given by 

(0|6.,(yjz;ie6l,(fci)at (fc)|0)(0|fo.,(fc3)a.,(fc)z;ie&L(^'e)|0) 
iVi = 2^ , =^ (63) 

\jvtp'e^kXkt (167r3(53(fc + fc3 -Pe)) (16^3J3(fc + fci -p'J) 

(for uie see Eq. (^), which after some algebra becomes 
Nx = (47ra) 



Pe A:;^/ \p'p^ k 



3 / Vi^e "-1 



2fc^ A:l(se) p^'(Se)W2fc' pl{s,) klise) 



yk+ fc+ p'/ J \ k+ p+ fc+ 

where we have defined a new object, 



(64) 



p'{s) = p^ + iseyp^. (65) 

Notation: i = 1, 2 only, s — ±1 only, s — — s, ei2 = —£21 = 1 and en — £22 = 0. 

We now discuss 9li and then simplify SMf^^ further. Recall Eqs. ( p^ and (p9[). We see that after the 
coordinate change defined by Eq. (pl|), in the rrip — > 00 limit, 



= eiX-T,)9{~X-^ + 0{a'))erX-^ + 0{a')] , (66) 



Ti = |k| + V(P - k)2 + m2 - me , (67) 
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where we have used the fact that the wave functions restrict |p| ~ ma. RecaU that we are always 
assuming ma^ ^ A <C ma. Thus, 9li can be simphfied: 



7L1 



^(A-Ti) . (68) 
From the form of Eq. (^), we see that this constrains the photon momentum to satisfy 

|k| < A, (69) 

to leading order in a. 

Note that the constraints coming from 9[^i, summarized by Eq. (|69|), require the photon momenta in 
SMl^ of Eq. (||) to satisfy 

k < pe,p'^ . (70) 

Thus, Eq. (|6^ can be simplified further. 



N' 



Ik Jpc J p'^ 

X (l)N' {p'^,V - p'J (f>*N, ipe.V ~Pe)(l)N {pcV - Pe) JJ^^^ 

In the TTip — > oo limit, — *■ Wp, and DENi becomes 



(71) 



(fe3=P.,fci=P^,|k|<A) 



DENi = 2mp(/3„-/3„, - |k|)(l + 0(l/mp)) , (72) 

where we have used + fc+ = 2|k|, valid for an on-mass-shell photon (all particles in a Hamiltonian 
approach are on-mass-shell). — /3„ is the binding energy of non-relativistic hydrogen defined in Eq. (^^, 
with numerical value Ryd/n^ for the bound-states. 

In the region of integration, |k| < A = ma^/a <C |p|, so after the coordinate change of Eq. ( pi] ) (recall 
Eq. (p3|)) we have 

iVi 4fc^^ 4fc-L^ 4 • p'^ 



47ra fc+TO TO^ 



^■(p-+p'--^-^)+0 (a^V^) . (73) 
K+m V mm/ 
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The rest of the integrand is even under k-^ ^k-^, so these terms in the last hne, odd in k-^, do not 
contribute. 

Putting it aU together, recalhng Eq. js^), we have 
SMIj^ ^ a f (pk 

N' 



- ^ » 4^1:7 W " - I"!) / ''^ I "V*™ (pO (P') 



fc+^ + fe'+m 



where we recaUcd Eq. (|4^), the relation between and ^^r. This is infrared ^ 0) divergent, but 
we must add diagram L2 of Fig. 2 to get the total low-energy shift. 

As previously mentioned. Diagram L2 of Fig. 2 arises from the sum of an effective second-order electron 
self-energy interaction and a counterterm defined such that the electron self-energy runs coherently. The 
result of this interaction is to add the following to the binding 



N 



X 0^, (p) (P) ^ ^ = . (75) 



Given the constraint |k| < A <C |p|, this becomes 

X (P) (P) ^ . (76) 

This is the famous subtraction that Bethe performed in 1947 |l| . In our approach it arose as a consequence 
of coupling coherence. 

5Bl2 is infrared divergent (fc+ — > 0) as is 5Bli. This divergence arises from the first two terms of 
A^i (the ones independent of p and p'). Noting that 

we have 

« fl 11 I (-k+^n) k^^ 

Pn - Pn' - |k| > , 
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and these infrared divergent contributions from the first two terms of A'^i cancel, leaving an infrared finite 
shift, 



X 0^, (p) 4)N (p) 2 



1 



1 



(3n - f3n' - |k| |k| 



lU (77) 



3 / 47r2 



E V M ^ " |k|) / d^P / d'p'^N (P') 0A^' (pO 



X (p) (/.AT (p) 



P P 



1 



1 



m2 V/3n-/5„' - |k| |k|. 
This last step followed after averaging over directions as dictated by rotational invariance. 
Eq. (ra) is easy to integrate, and our final result for the low-energy photon contribution is 



6Bt 



2a 
Svrm- 

2a 



iPn' - f3n) log 



N' 



X + f3n'- Pn 



■5]'(/3„'-/3„)log 



A" 



I3n' - A 

A 



|(0Af|p|(?!>7V')r 



|(0Ar|p|0JV')r 



(78) 



(79) 
(80) 



where in this last step we recalled A ^ raa^ . Note the A-dependence in the result. This will cancel after 
we correctly add the contributions coming from high energy intermediate photons, which now follows. 

B. High energy contribution 



The high energy shift arises from three sources which are shown in Fig. 4. These are first-order BSPT 
shifts due to third and fourth order effective interactions (see Appendix B) . The net result of these three 
diagrams is 



27r2q2 27r2q2 
where q is the exchanged momentum of the electron, and 



{1 + 5Vh) , 



(81) 



6Vh = 6Vhi + 6Vh2 + 5Vh3 



(82) 



with 



5Vhi = 



k+ 6 



{pt - k+) Nhi Ohi 
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{V+f {V+f 



(M2 - M2) (M2 _ (Mf - M2) (Af^2 _ ^^,2^ 



(83) 



5Vh2 = 2 



/f3 — 



The factors ±5 in front arise from the form of the Bloch transformation (see Eq. (135) of Appendix B). 
The vertex factors are given by 

Nhi = (A^i)(,,^p,_fc,fc3^p^_fe) , (86) 
Nh2 = (A^i)(fcj^p/„fc^fc3^p^_fc,p^^p/j , (87) 

Nh3 = (A^l)(fei^p,_fe,fe3^p^_fc,p^^p^) ' 

where A^i was defined in Eq. (Bih. The free state masses are given by 



= Vp' + Jp2 + m2 , (89) 



K = sJp'' + nil + Jp'' + ml, (90) 



M = |k| + v/(p - k)2 + m2 + w'p2 + to2 ^ (91) 



M' = |k| + v/(p' - k)2 + m2 + ^p'2 + ml . (92) 

The Bloch transformation constrains the free masses of the states. As discussed before, the "L" 
restrictions in Fig. 4 can be removed given A 3> ma^. However, the "H" restrictions lead to important 
constraints given by the 9h factors above, which we now discuss. They constrain the free masses to 



satisfy (recall Eqs. (29)-(pa 



A < A/ - me - mp < A , (93) 

A < A'f - TOe - TOp < A , (94) 

where M and M' are defined in Eqs. (^ij) and ( p^ ) respectively. 

As already mentioned, for convenience of calculation, we will divide this high energy contribution into 
two regions, A < |k| < 6 and & < |k| < to (region one and region two respectively), with ma <C 6 ^ 
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m^/a ^ m. Recall, ma^ <C A ~ ma^/a ^ ma. We now show how this division into these two regions 



arises as a result of the constraints of Eqs. (|93| ) and (|94| ) 



In this first region, |k| <ti m, and Eqs. ( |93| ) and ( |94[ ) become 

A < |k| + ^ |k| < , (95) 

A < |k| + ^ |k| < 6 , (96) 

which is as we have already stated (recall that we always assume nip — > oo and drop the 1 / corrections 
since we are just after the dominant shift). 

The analysis of the second region is slightly more complicated because |k| ^ ma, and near the upper 
limit |k| ~ m. Since |k| 3> ma in this region, Eqs. (BS) and (p4) both become 



b < |k| + Vk2 + m2-TO < A. (97) 

This is just a linear constraint, 

b(^) < |k| < Ifii^l , (98) 

\2m + 2bJ - ''-2\^A + my ^ ' 

which, since we choose b <^ m, becomes 

'^""4 (It!:)- 

The electron self-energy renormalization is performed in this paper, but we do not deal with removing 
the full A-dependence. A full analysis of this dependence requires a complete 4th order calculation, which 
is beyond the scope of this paper. We cut off the photon momentum at the electron mass, and proceed. 
Note that from Eq. (p9[), this choice corresponds to A^ = 2m^. The point of calculating these high photon 
energy contributions is to show that our results are independent of the effective cutoff, A. 
Taking a sample denominator we have 

(M2 - m2) ^ {M, + M){Mo - M) « 2mp - |k| - ^^^^^ ~ - 2mp|k| , (100) 

in the first region; and 
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(M2 - M^) = {Mo + M){Mo - M) « 2mp (m - |k| - + m^) « -2mp (^|k| + (101) 
in the second region. 

Using these previous formulae, including — > rup as nip — > oo, Eqs. (^3|)-(p5|), after summing, 
become 

in the first region (the "prime" on 5Vh signifies the first region); and 

" 9^' r . n /""^Ikl A |k| |k| 



o ( § n > (103) 



in the second region (the "double prime" on SVu signifies the second region) . In the second region since 
the photon momentum is not necessarily smaller than the electron mass, we have kept two terms in 
the |k|/m expansion of the integrand. In the Odkl/m) term we have introduced two constants, c„ and 
Cd, which denote numerator and energy denominator corrections respectively. Hereafter we set c„ — 1 
and Q = 1, as given by the theory. Note that in combining 5Vhi, SVh2 and SVh3, many cancelations 
occur; most noteworthy, each contribution is individually infrared divergent (fc+ — > 0), but in the sum 
the divergences cancel. These final equations are easily integrated, and we have 

SVh = -^^logfl) , (104) 



27r TO^ \A 
zn 771^ V / 

In the second region note that the Odkl/Tn) terms coming from numerator and energy denominator 
corrections cancel, leaving the 0(1) piece alone. The combined high energy contribution is 

5Vh = SVi, + 5Vii^-^^ log f"^) , (106) 



27r 77l2 "° V A , 

which is independent of b, as required for consistency. Recall that q = p' — p: the difference between the 
final and initial electron momenta. 
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From the definition of 6Vh (see Eq. (81)), we see that this correction shifts the energy levels an amount 



Combining this with Eq. (|l06[) gives 



2 ; <^A' (P) 



(107) 



SB 



H 



log (^) / d^p d^p'r. (P') f ^^pl^y ) 0. (P) 



„2 



log(^^) d3p0Ar(p) 



(108) 

(109) 



where in this last step we averaged over directions and noted that the wave function at the origin is real. 
For more details on this averaging over directions see Appendix C. 

C. Total contribution 



In this section we combine the results of the last two sections for the low and high photon energy 
contributions, and perform the required sums/integrations to calculate the total shift between the 2Si 
and 2Pi energy levels of hydrogen. 

Adding Eqs. (|^) and ( [l09| ) gives for the total shift 



SB = SBl + SBh = 



2a 



■5]'(/3„' -/?n)log 



N' 



A 



|('/>Ar|p|0Af )r 



a" I m 

log , 

A 



d^p (j)M (p) 



For the second term we have 



The (27r)^ factor arose because of our normalization choice (see Eq. (p|)). 



(110) 



(^j d^p<P^{Y>)j = ((27r)i0^.(x = O))' = M^(^)'<S,o. (Ill) 



The first term of Eq. (110) is the famous Bethe log and must be calculated numerically, summing over 
all bound and continuum states. Following standard convention we define an average excitation energy, 
^(n,0. 
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^Og[^^]J2^{Pu'-f3n)\{cbN\p\ 



Ryd 



if^n' - Pn) log 



N' 



N' 



Ryd 



\{^n\p\M\ ■ (112) 



The sum on the left is evaluated by standard techniques {He — /{2m) — a/r ), 



5^'(/3«'-/3n)|(0jv|p|0. 

1 , 



2 1 

'N']\ = Tii'l^N \[p,Hc]-p + p- [Hc,iy]\(j)N) 



N 



ON 



5W 



p-, — iV {—a/r) 



\ (^AT |(-*)2V2 {-a/r)\<t>N) - -\{-i)\-a){-^'K) {^n \S' (r)| <Pn) 



2a (^YS^ 



\ n / 



1,0 



(113) 



This vanishes for I ^ 0, but the average excitation energy, /3(n, I), is defined (it is just a way to catalogue 



the numerical sum on the right of Eq. ( 112 ), the quantity we need to know) with the sum on the left 
hand side set to its value for 1 = 0. In summary, P{n,l) for all states is defined by 



Ryd 



i((/)jvipi0iv')r 



(114) 



Without further ado, this sum has been evaluated by R. W. Huff M. His results for the ?i = 2 levels are 



/3(2,0) = 16.63934203(1) , 
P{2, 1) = 0.9704293186(3) Ryd , 



(115) 
(116) 



where the figures in parentheses give the number of units of estimated error in the last decimal place (R. 
W. Huff's estimates). 

Combining the results: 



5B2S, 



6B' 



2a 
37rm^ 

a^Ryd 
37r 
2a 



log 



I ^12 (^"^Y 



a^ 1 ( ™^ (2"')'^ ^ma\^ 
GTT^m^ ^ V A 



3 ,^„,s3 

TT V n / 



log 



/3(2,0) 



2Pi 



\ n / 



37rm2 '^V/3(2,l) 

a^Ryd , / Ryd 
■ log ^ 



3^ "°V/3(2,1) 
Note the cancelation of the A-dependence. Thus, the Lamb shift is 



(117) 



(118) 
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^B,^^, = m^s. - 5B2P, = — ^ log ^.i. ' 119 
5 3 dvr \RydP{2,0)/ 

= (1047 - 4) MHz {2T:h) = 1043 MHz {2T:h) , (120) 



where we use Ref. 1^ and the average excitation energies of Eqs. ( |ll5|) and (116). Note that the 2Pi 
shift is only about one half of a percent of the 2Si shift. 

V. SUMMARY AND DISCUSSION 

In a light-front Hamiltonian approach, we have shown how to do a consistent Lamb shift calculation 
for the n = 2, j = l/2 levels of hydrogen over the photon energy scales 

<-> ma^ ^ A ^ ma ^ b ^ m , 

with the choices ma^ ^ A ^ ma and ma <C <C m. In a consistent set of diagrams we showed how A- 
and 6-dependence cancel leaving the dominant part of the Lamb shift, 1043 MHz. For completeness, the 
n = 2 spectrum of hydrogen is shown in Fig. 5. 

If we compare the three regions we see the following results (we only compare for the 25*1 shift since 
the 2Pi shift is negligible within our errors): 

(0<|k|<A) ^^....-^log(^^:^)-46MHz(2..)^4%, (121) 
(A < |k| < b) SB^^^, - log (i) - 667 MHz (2^;i) ~ 64% , (122) 



3^ VA. 

(b < |k| < m) SB^ , - log f ^"j - 334 MHz {2Trh) - 32% , (123) 

where we used A = may/a and b = m^/a, consistent choices used throughout this paper. As expected on 
physical grounds (see the Introduction), photons with momentum 

|k|-l/ao-ma, (124) 

couple the strongest to the hydrogen system. As seen above, the effects of photons of this momentum 
amounted to about 2/3 of the Lamb shift, the dominant part of this experimentally observed shift. 

In this paper, the one loop electron self-energy renormalization was performed. The complete one loop 
renormalization was not needed to obtain the dominant part of the Lamb shift. Our answer, 1043 MHz, 
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turned out to be accurate. However, to obtain more precision, the full one loop renormalization must be 
performed of course. Also, each of our five diagrams (of Figures 2 and 4) were infrared {k^h^ton — * 0) 
divergent. However, both the sum of the two low photon energy diagrams and the sum of the three high 
photon energy diagrams were infrared finite. 

The state of the art for the bound-state problem in a light-front Hamiltonian gauge theory in four 
dimensions has been advanced in this paper. In applying these methods to QCD, the general computa- 
tional strategy described in the Introduction docs not change. However, gluons carry color charge and 
interact strongly at low energies, thus the answers to the three questions posed in the initial paragraph 
of the Introduction change drastically. For a constituent picture to emerge the massless gluons must be 
confined, so that it costs energy to add a low momentum gluon to the system.]^ It has been shown that 
the second order effective interactions (including the very important first order instantaneous-gluon po- 
tential) are confining 1^, which is promising. Given confinement, we can lower the effective cutoff below 
the gluon production threshold perturbatively and obtain a constituent approximation. As in QED, we 
can not lower the effective cutoff below the non-perturbative bound-state energy scale. Thus Eq. (^) of 
hydrogen in QCD becomes 

AqCD ^ a ^ Egiuon ^ Mgiueball /2 . (125) 

Since Aqcd ranges from 200-400 MeV (depending on the renormalization scheme that is chosen) and 
Mgiueball ranges from 1500-1700 MeV, this constraint can be satisfied and it becomes plausible to attack 
QCD by the same computational strategy that was outlined and carried out for the Lamb shift in QED 
in this paper. 



The energy of this confined low momentum gluon can be interpreted as an effective gluon mass if it is convenient. 
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APPENDIX A: LIGHT FRONT CONVENTIONS 



In this Appendix we write our light-front conventions for the electron, proton, and photon system. 
• V^ = V°±V^ where is any four vector. 













-2i 








; 7 = 






2i 











e+ 



1 



cr' 

a* 

1 



i = 1,2 ; CT* are SU{2) Pauli matrices. 





'4)± = A±v ; ip = + '4' 



; A_ = ^7°7- = 




1 







e for electron, p for proton. 





^ [{a' {id' + eA')+ime)^e] 




J+ = 2e {i% - el^e) ; e > 
A+ =0 
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In momentum space the field operators are expanded as (at = 0): 



A' 



s=±l 



with 



Xi = 



Jq 


f h.c.) , 


Xs V 




+ 4(we ) . 


Jp 








Xs V 








Jp 








-1 








— ((5,;,i + i di^2) 

v2 






(^) 

























Note that s = —s. Also, we are using the shorthand 

The fermion helicity can only take on the values ±1/2, however we define: hj, = s/2; therefore, "s" can 
only take on the values ±1. The commutation (anti-commutation) relations and free Fock states are 
given by: 

[ax{q),a\,{q')] = Wir^+S^q - q')S^y , ( 6\p) = 5\p^)5{p+) ) , 
{bs{p),bl,{p')} = {4(P),4(P')} = 167rV'5^(p-p')'5..' , 
{B,{p),Bl{p')} = {D,{p),Dl,{p')} = lQi,^p+5%p~p')5ss' , 

{piSi\p2S2) = 167T^p^6^{pi ~P2)Ss^s2 , \PlSl) = &Ii(Pl)|0) , ctc . 

APPENDIX B: BLOCH TRANSFORMATION 

In this Appendix we discuss a derivation of our effective Hamiltonian via a Bloch transformation Q . 
We use the Bloch transformation to separate the low and high energy scales of the problem and derive 
an effective Hamiltonian acting in the low energy space alone with an identical low energy spectrum to 
the bare Hamiltonian. In this Appendix, we closely follow §IV of Ref. |7j, where a discussion, including 
the original references, and derivation of a general effective Bloch Hamiltonian can be found. 
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We start with a bare time-independent Schrodinger equation: 

Ha\^^) =E\^^) . (126) 
Then projection operators onto the low and high energy spaces, Pl and Ph respectively, are defined, 

/\2 I -p±2 \ 

Pl = 0{^T 'M' (127) 

Pl + Ph = ^(^T '^l ' (129) 

where 9{x) is a step function. Then an effective Hamiltonian acting in the low energy space alone with an 
equivalent low energy spectrum to H\ is sought. A and A are the bare and effective cutoffs respectively 
with A < A.^ V = (V^-.V'^) is the total momentum of the hydrogen state, h is the free Hamiltonian of 
the hydrogen system of Eq. (||). 

Proceeding, a new operator, TZ, is defined that connects the Pl and Ph spaces: 

Ph|*a) =7^PL|*A) . (130) 

More explicitly {J2n = !)• 

n,m \ I I / 

which shows that TZPhI'^a) ~ and 7?.'''Pl|^a) = 0, etc. For the construction of 7?., see below Eq. (4.4) 
in Ref. §. 

This leads to the following time-independent Schrodinger equation for the effective Hamiltonian 

Hxl'i'x) ^E\<i>x) . (132) 

E is the same eigenvalue as in Eq. ( |126D . The state \^x) is a projection onto the low energy space with 
the same norm as Iv^a) in Eq- ( |l26|) : 



°The shorthands — > A and — > A are often used when it does not lead to confusion. 
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i$a) - yrT^PLi*A>; (133) 



H\ is a hcrmitian effective Hamiltonian given by 



Hx = , ^ , {Pl + Ha {Pl + n) ^ ^ , . (134) 



Note that Hx acts in the low energy space alone. To summarize, Hx of Eq. (134) is guaranteed to have 
the same low energy spectrum as the bare Hamiltonian, Ha] also, after diagonalizing Hx, the bare state, 
I^-a), of Eq. ( p^ , if desired, is obtained through Eqs.(|l3|) and (|l3^).[| 



Defining Ha ~ h+VA, where h is the free field theoretic Hamiltonian and va are the bare interactions,^ 
to third order in va, the effective Hamiltonian is given by 

/ irr |,\ / , iLX , 1 {a\vA\i){i\vA\b) {a\vA\i){i\vA\b) 
{a\Hx\b) = {a\h + VA\b) + - I ^- + ^- 

^ \ ai hi 

{a\vA\i){i\vA\i){3\vK\b) , {a\v A\i) A\j) {i\v A\b) 



2 ^ V AaiAaj AbiAb,- 



( {a\vA\i){i\vK\c){c\vA\b) , {a\v a\c) {c\v A\i) {i\v aW \ , . 



o{^{vA)V{eMgh-eio..f) 



A-ia = Si — Sa, with h\i) — ei\i). We are using |a), • • • to denote low energy states (states in Pl) 
and • • • to denote high energy states (states in Ph)- See the already mentioned Ref. for a 

description of an arbitrary order (in perturbation theory) effective Hamiltonian and also for a convenient 
diagrammatic representation of the same. 



'^This state, |^a)i will span the whole space, but will correspond to the respective low energy eigenvalue. The bare 

states that correspond to the respective high energy eigenvalues can not be obtained from the effective Hamiltonian, Hx\ 

we must of course use the bare Hamiltonian, Ha, to accomplish this task. 

® /i is written in terms of renormalized parameters, and it is convenient to define va = v+5v/^ , where v is the canonical field 

theoretic interactions written in terms of renormalized parameters and Sva are the counterterms that must be determined 

through the process of renormalization. See Eqs. and for the canonical Hamiltonian of the hydrogen system. 
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APPENDIX C: AVERAGING OVER DIRECTIONS 



In this Appendix we calculate the following Coulomb matrix element: 



/± = 1 d'p J d^p' ^^(p')l^i-^J^0Ar(p) , (136) 



to verify the step taken from Eq ( |108| ) to (|109| ) in the paper. It is useful to define another integral, 



d^p / dV 0^(p')^r^-^0iv(p) • (137) 

(p-p) 



Now note that 



/ = + I d^p f d^p' 0^(p')0Ar(p) = ^ {^y^''° ' (138) 



where in this last step we recalled Eq. (Ill) and the fact that the wave function at the origin is real. 
For I = 0, the wave function satisfies 

(t^nfifiip) = (t>nfi,o{\P\) ■ (139) 

Thus, by symmetry, for I — 0, 

/.4/=^(^)\ (140) 
3 37r V 71 / 

which from Eq. (|l38| ), for 1 = 0, gives 



3 37r \ n / ' 



(141) 



For I ^ 0, first note that 7 = 0. Thus, for I ^ 0, 

/± = -h ; (142) 
we will calculate Iz below which then implies Ij^ . Next note that in position space 

/ = ^2^^ j d''x4>*^{^)(v^l-^cj,N{^) , (143) 
using V^i^ = —At: 5^ (x). Thus, for Z 7^ 0, in position space 
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h = -27r2y d?x 4>*N (x) (^'^j (x) . (144) 

Note that there is no |x| ambiguity in this previous equation because for I ^ 0, the wave function 
vanishes at the origin. Carrying out the derivative gives 

/. = -27r2 / d?x 4>*M (x) ( "l + ^^'/N" ) (x) . (145) 



This matrix element was performed in the first appendix of Bethe and Salpeter's textbook ||lo|] . We use 
two of their formulas, (3.26) and (A.29).0 Eq. (|l45| ) integrated gives 

h = -2tt^V=^ c{l,mi) , (146) 

with 

— =■ 1 / ma \ 3 



- 



Z(/ + l)(/ + i) V 



(147) 



- + (2/ + 3)(2/-l) j • 
Thus, recaUing Eq. ( |142D , our result for ^ ^ is 

= 27r2 ^ c(Z,m/) . (149) 

For Z = 1, /j^ is not zero, so what is going on? The answer lies in the fact that we really want to 
take matrix elements in the \j,mj) basis not the \mi,ms) basis, and based on rotational invariance, our 
results should be independent of rrij. To proceed, note that the interactions we considered in this paper 
conserved nis and our matrix elements were independent of mg. Next, note that the result, Eq. (149), is 
even under m; — > —mi. Given this, the Clebsch-Gordan coefficients for the 2Pi states imply 

1 ' 

(j = l/2,m,|y|j = l/2,m,) = {mi\V\mi) , (150) 

mi=-l 



where 1 = 1. Now note that Ij_ given by Eq. (149) averaged over mi vanishes 



^ E ^^=0, (151) 



21 , . 

mi= — l 



^Warning to the reader: In this text, they use atomic units, h = c = m = a = 1, so m and a have to be placed back 
into the formulas. 
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where we used 



21 



(152) 



nii = -l 



an obvious result after the answer is known. This result (Eq.(151)) was used in the step that led from 
Eq. (|108|) to Eq. (|109|) in the paper, and this Appendix is now complete. 
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Figure captions 

Figure 1: The effective interactions that add to give the Coulomb potential. "H" implies that the 
photon energy is greater than A. "L" implies that the electron kinetic energy is less than A. We choose 
ma'^ ^ A <C ma; these "H" and "L" constraints can thus be removed to leading order. 

Figure 2: The low energy contribution of §IV.A. Diagram LI represents the shift arising from treating 
photon emission below the cutoff A in second-order BSPT, where the intermediate electron-proton are 
bound by the Coulomb potential. Diagram L2 is an effective self-energy interaction (plus counterterm) , 
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arising from the removal of photon emission above the cutoff A, treated in first-order BSPT. /3(2, 1) is the 



average excitation energy of the n = 2 state; see Eqs. ( |f 15| ) and (|f f 6| ) and the discussion above them for 
details. 

Figure 3: The sum of an effective self-energy interaction arising from the removal of photon emission 
above the cutoff A and a counterterm. The counterterm is fixed by coupling coherence. The result is the 
interaction in diagram L2 of Fig. 2. 

Figure 4: The high energy contribution of §IV.B. These are third and fourth order effective interactions 
treated in first-order BSPT. These effective interactions arise from the removal of photon emission above 
the cutoff A. '6' is an arbitrary scale, required to satisfy ma <^ b <ti m, that was introduced to simplify 
the calculation. Note the 6-independence of the result. The total contribution is a sum of Fig. 2 and 
Fig. 4. Note the A-independence of the combined result. 

Figure 5: The n = 2 hydrogen spectrum: fine-structure, Lamb-shift and hyperfine structure. F = 

L + Se + S„. 
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SELF-ENERGY RENORMALIZATION 
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THE HYDROGEN SPECTRUM: n = 2 
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